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Abstract. - We study the steady state of a multiply-connected system that is driven out of 
equilibrium by a sparse perturbation. The prototype example is an iV-site ring coupled to a 
thermal bath, driven by a stationary source that induces transitions with log-wide distributed 
rates. An induced current arises, which is controlled by the strength of the driving, and an 
associated topological term appears in the expression for the energy absorption rate. Due to the 
sparsity, the crossover from linear response to saturation is mediated by an intermediate regime, 
where the current is exponentially small in y/~N , which is related to the work of Sinai on "random 
walk in a random environment" . 



The transport in a chain due to non-symmetric transi- 
tion probabilities is a fundamental problem in statistical 
mechanics [TH5]. It can be regarded as a random walk in 
a random environment. The seminal observation is due to 
Sinai 0: considering a chain of length N, the random- 
ness implies a buildup of an exponentially large potential 
barrier exp(\/iV), and consequently an exponential sup- 
pression of the current, reflecting a sub-diffusive [log(t)] 4 
spreading in time. 

We would like to explore how this picture is modified if 
the chain is replaced by a configuration with a non-trivial 
topology, such as a ring, accounting for: (1) unavoidable 
telescopic correlations between the transition probabili- 
ties; (2) sparsity due to log-wide distribution of the tran- 
sition rates as in glassy systems. 

Let us expand on the roles played by these two ingredi- 
ents: In the physical model that we would like to consider 
(Fig.[T]) the transition rates are not independent random 
variables. Rather they arc determined by differences of un- 
corrected on-site energies. These differences are strongly 
correlated: their sum does not grow with N. Consequently 
the theory of Sinai does not apply. We shall see that in or- 
der to witness a "Sinai regime" we have to introduce into 
the model an additional ingredient: one may call it glassi- 
ness or sparsity. These terms emphasize complementary 
characteristics of the assumed log- wide distribution of the 
transition rates: the distribution span several decades; ac- 
cordingly the vanishingly small values constitute the ma- 
jority; while the large values constitute a minority; which 



is reflected by having a median that is much smaller com- 
pared with the algebraic average. 

Our focus is on the non-equilibrium steady state (NESS) 
global current / that circulates the whole ring, and on the 
associated energy absorption rate (EAR). Let us present 
some numerical results that clarify the physical picture 
and motivate the subsequent analysis. We consider a ring 
that is composed of N sites (Fig^TJ). The ring is weakly 
coupled to a bath that has temperature Tb- In the ab- 
sence of driving the average current I is zero. The driving 
is modeled as a stationary noise source that has an in- 
tensity e 2 . The driving breaks detailed balance, leading 
to a non-vanishing affinity along the ring. This affinity, 
which we call stochastic motive force (SMF), labeled as 
£q, induces a circulating steady-state current, see Figj2l 
For weak driving one observes a linear response behavior 
/ oc E(5 a e 2 . For very strong driving / saturates. 

Both the linear response for weak driving and the sat- 
uration for strong driving have a relatively simple expla- 
nation (see details later). But the crossover from linear 
response to saturation is more subtle, because it is related 
to the distribution of the transitions rates. If the distribu- 
tion is log-wide we call the system sparse, implying that 
most of elements are much smaller compared with the av- 
erage. Due to the sparsity one observes that there is an 
intermediate region where &q ~ VN , while the current 
becomes exponentially small, exhibiting fluctuations as a 
function of e 2 . 
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be regarded as the entropy change of the bath during the 
(m — * n) transition [TJ. Consequently it is natural to de- 
fine the "potential variation" between two points x\ and 
X2 along a line segment as follows: 



£(xi x 2 ) 



£{x)dx 



(4) 



Given a loop, one defines the SMF (or mesoscopic affin- 
ity [S|) as follows: 



£q = y^£(x) — (p £{x)dx 



(5) 



Fig. 1: A ring made up of N isolated sites with on site energies 
e„. The ring is coupled to a heat reservoir (represented by 
the blue "environment") and subjected to a noisy driving field 
(represented by the red circle) that induces a current in the 
ring. In the numerical tests the energies occupy a band of 
width A = 1, and the temperature of the bath is Tb = 2. The 
driving source induced rates are log-box distributed over 
8 decades, while the bath induced rates are all with w^- — 1. 



Outline.— We introduce the Ring model and clarify 
its relation to the standard paradigm of NESS analysis. 
We derive expressions for the SMF, for the current, and 
for the EAR, illuminating their dependence on the driv- 
ing intensity. With regard to the EAR, we highlight the 
manifestation of a topological term. 

Sparse networks.— Consider a general rate equation 
for the occupation probabilities p n , namely 



dpn 
dt 



J nrnJr / m 



(i) 



where w nm is the transition rate from m to n. We can 
regard it as describing a network that consists of "sites" 
connected by " bonds" . Specifically we consider later a ring 
that consists of N sites (FigfTJ). With regard to the bond 
x = (m n), that connects site m to site n, we define the 
coupling w(x) and the field £{x) as follows: 



£(x) 
w(x) 



In 



J nm "-'mn I 



,1/2 



(2) 
(3) 



We say that the system is sparse or glassy if either w(x) or 
£(x) of the connecting bonds have a log- wide distribution. 
This means that there is a small fraction of bonds where 
the coupling or the field is strong, while in the overwhelm- 
ing majority it is very small. 

We can regard £(x) as the "potential difference" across 
the bond x = (m~» n). The reason for this terminol- 
ogy becomes obvious if one considers bath induced tran- 
sition [see later discussion regarding Eq. (|7J]: Then it 
equals — (E n — E m )/Tg- Note that optionally £(x) can 



The summation above is over the bonds x along the loop, 
which becomes an integral in the continuum limit. For a 
detailed balanced system the SMF is zero for any loop. 
Otherwise the system relaxes to a NESS. 

For the subsequent analysis we define "the effective po- 
tential barrier" along a segment as the maximal potential 
variation (see how it looks in Fig|3|): 



£ n — maximum | \£(x± x%] 



(6) 



Referring to a Ring, it is important to realize that <?n 
cannot be smaller than £q. If the £{x) were totally mi- 
correlated both the maximal potential variation and the 
SMF would be proportional to y/~N. 

NESS paradigm.— In the physical problem the net- 
work consists of N sites, with on-site energies E n . The 
transition rates w nm from site n to site m are induced by 
a driving source that has an intensity e 2 , and by a bath 
that has a temperature Tb- Namely, 



it', 



(7) 



where w e nm = 
balanced with 



while the bath 
w n m /wLn = exp[-(E n -E n 
can define a symmetrized rate as Wnfn - 
Then it follows that 



(w, 



is detailed- 
)/T B }. One 



,<3 - 



2<m 



I + exp 



(8) 



Our setting as described above is formally a special case 
of the common non-equilibrium paradigm for a system 
that is coupled to two heat baths. The driving source 
is like a bath that has temperature Ta = oo, while the en- 
vironment is a bath that has a finite temperature Tb < oo. 
As discussed in Ref. [1] the generalization of the detailed 
balance condition requires 



{Qa, Qb) 







exp 



Qa 



(9) 



Where Qa and Q b are the the amounts of energies that are 
transferred from the baths to the system. In our setting 
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the two baths arc independent of each other and therefore 
any event is cither with Qa = E n — E m and Qb = 0, or 
with Q B — E n — E m and Qa = 0. 

Microscopic temperature.— In the absence of driv- 
ing the transitions that are induced by the bath sat- 
isfy detailed balance, and accordingly the steady state 
of the system is canonical with occupation probabilities 
p n oc exp[— Eu/Tb]- Once we add the driving this is no 
longer true. Still, there is a well defined NESS, so we can 
formally define a microscopic temperature for each tran- 
sition separately via the formula 



Pn 

— = exp 

Prn 



T 



(10) 



Unlike a canonical state, here we may have a wide dis- 
tribution of microscopic temperatures. Furthermore, in a 
NESS the local temperature can be negative, i.e. the occu- 
pation of a higher level can be larger than the occupation 
of a lower one. 

The Ring model.— As a specific example we con- 
sider a ring with random on-site energies E n £ [0,A], 
and near neighbor transitions. We use the notations 
A„ = E n - £7 n _i, and Wjf = uy„_i, and = w„_i,„, 
and Wn = (u)jf +wfc)/2. The superscripts /3 and e are 
used in order to distinguish the bath and driving source 
contributions. It is useful to notice that in the high tem- 
perature limit Eq.© takes the form 



it 



A, 



2T B 



(11) 



The rate equation Eq.([T]) is formally identical to an elec- 
trical network where the p n are like voltages, and the w nm 
represent conductors. Inspired by this formal analogy, we 
regard the conductivity of the network as an "average" 
transition rate, which we call w. Consequently for a net- 
work with near-neighbor transitions, which is like having 
"connectors in series" , we get 



N E 



w(x) 



(12) 



with similar definitions for w° and w e . It is now natural to 
define a dimensionless driving intensity e 2 and dimension- 
less coupling parameters g n , such that the latter reflect 
the relative exposure of the bonds to the driving: 



w 
K 



9n£ 



[definition of e 2 ] 
[definition of g n 



(13) 
(14) 



Sparsity means that the couplings to the driving source 
are log-wide distributed. In the numerics we have as- 
sumed log-box distribution of the g n over several decades. 



Namely, we the values \n(g n ) were generated such that 
they form a uniform distribution in the range [g m in, flWx]- 
As for the bath: the were assumed identical. Accord- 



ingly note that l/g n 
ance is Var(</„) = g n 



fJ 



1, while 1 < g n <C g% 

2 



The vari- 



Estimating the SMF.— In the presence of driving 
the SMF is non zero: 



Sr, = In 



where Y[° is the product of all the N anticlockwise rates, 
and Y[ is similarly defined. If Tg 3> A we get 



[O] 



(15) 



£n = In 



n„«+^) 



E ln 

n 

E 



i - 



1T B 



1 



l + g n e 2 



1 + Ifej 

An 

T B 



(16) 



Recall that we have A n = 0. Additionally we define 



A (o) 



A (°°) 



E Sn^r, 



„ 9n 



± 



± 



2iVVar( ff ) 



1/2 



A 



1/2 



(17) 



2iWar( ff " 1 ) A (18) 



The RMS-based estimate of the sums follows from 
the observation that, say, A(°) can be rearranged as 
^2 n {g n +i — gn)E n , which is a sum of N independent ran- 
dom variables. Consequently we get for the SMF the fol- 
lowing approximation 



J_ _ A (~)/ e 2 , 
-<[±]A«, 



e l < l/g, 



e > l/.g min 
otherwise 



(19) 



where AM = iV 1 / 2 A, and the ~ implies that the the result 
exhibit fluctuations as e 2 is varied. Looking in FigJ^l at 
the plot of we notice that there are dips that indicate 
that the SMF changes sign. These sign reversals depend 
on the specific realization of g n and A„. 

The NESS current.— The rate equation for nearest 
neighbor transitions is 



(20) 



This set of equations is redundant because of conservation 
of probability. At steady state p n = 0, and there is some 
current I in the system. So we can write the equivalent 
non- redundant set of N+l equations 



WjfPn-l ~ W^p n = I, 



Pn = 1 



(21) 
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Fig. 2: We consider the Ring of Figfflwith iV=10 3 sites. The 
absolute values of the SMF (solid curve) and the current (thick 
dashed curve) are plotted as a function of the scaled driving 
intensity. The vertical dashed lines l/g ma x and l/p m in bound 
the Sinai regime. Due to the sparsity (g m in <C gmax) it ex- 
tends over several decades. The solid horizontal line is the 
estimates for the SMF Eo. (|19[) . The two straight dotted lines 
(labeled as Eq.B) are the estimates for I based on Eq. (|25ll . 
while the dashed line (labeled as Eq.A) is the global approx- 
imation Eq. (|24[) . The diamonds indicate valued for which we 
plot the potential landscape in FigO 



where the denominator has a sum over bonds x and target 
sites n. In each term the bond x is disconnected to open 
the ring. The product rj^^") denotes the multiplication 
of the N— 1 rates that are directed from the cut point 
towards the target site n. Note that the TV— 1 rates belong 
to the two paths that lead from x to n. In continuum limit 
notations the expression can be written as 



3 £o/2 _ p-fo/2 



V i fdx'e £ (^')/ 2 

t-~/x w(x) J 

Which can be roughly estimated as 



(23) 



/ ~ — w exp 
N y 



2 



(So 



(24) 



This rough estimate is tested in Fig[2J and is in fact quite 
satisfactory. 

Let us discuss in more detail the current for very weak 
and very strong driving. In both limits the SMF becomes 
very small, the Sinai factor cxp[] becomes of order unity, 
and the sinh() can be approximated by a linear function. 
Accordingly we get 



1 

N 



-[A(°)/T B 
[A(°°)/Tb]w 



linear regime 
Saturation 



(25) 
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Fig. 3: The Ring of Fig[2]is considered. The potential difference 
£ (0~> x) is plotted against x, for two values of the driving in- 
tensity. The vertical lines correspond to the maximal potential 
variation £ n . The SMF values £q = £ (0 ~» N) are indicated 
by diamonds. 



This set of equations can be solved for the current using 
elementary algebra, or alternatively using the network for- 
malism for stochastic systems [SHlOj. The result may be 
written in compact notation as follows: 



t n 



[O] ~ [O] 



(22) 



As evident from these expressions, and as implied by the 
numerics, the direction of the current can change as the 
strength of the driving is increased, hence the dips in |7| 
in Fig|4] The small e result is independent of , in spite 
of the bath dominance. The fingerprints of the bath show 
up only for strong driving, where the result becomes sat- 
urated, independent of w e . The saturation by itself could 
have been anticipated: it is due to having SMF oc 1/e 2 
with rates oc e 2 . 

In the intermediate regime \£ \ - N 1 / 2 . The maximal 
potential variation £ n is of the same order of magnitude, 
but always lager, typically by some factor of order unity. 
Hence the current becomes exponentially small in \fN, as 
in the model by Sinai. Sparsity is the crucial requirement 
in order to observe this intermediate Sinai regime, other- 
wise there is merely a crossover from the linear response 
regime to the saturation regime. 

The EAR formula.— By construction it is obvious 
that for zero driving intensity the system is detailed bal- 
anced with the temperature that is dictated by the bath: 
hence Eq. (fit))) is satisfied with T nm = Tb ■ Also for non 
zero driving the system is detailed balanced, provided the 
current is zero. This is implied by Ea. (|21j) . Consequently, 
in the absence of current, Ea. (fT0| is still satisfied, but the 
microscopic temperature of the nth transition is 



1 

-L n 



l 



In 



"'1 
w-„ 



(26) 
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For Tb 3> A one easily obtains the following practical ap- 
proximation: 



1 




1 
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T (o) ~ 


1% 




.l+5n£ 2 . 



1 

T~B~ 



(27) 



One observes that in the absence of current the system is 
locally heated, with microscopic temperatures T^ > Tb 
that are non-uniform due to the dispersion in the cou- 
plings. 

Let us see what happens if the current is non zero. 
Again we can try to use Ea.(|10p in order to define a set 
of microscopic temperatures T n . One should notice that 
unlike T ( n\ some T n might be smaller than Tb or even 
negative reflecting a probability occupation inversion (in 
the sense of Laser physics). For the following analysis it 
is useful to observe that Eq.(|21|l provides an explicit ex- 
pression for the occupation difference: 



Pn-l —Pn 



2p n 



w^+w-, 



(28) 



where p n = (p n -i +Pn)/2 is defined as the average oc- 
cupation of the sites in the nth bond. Note also that 
w<n—Wjt originates exclusively from the bath, because the 
driving induces symmetric transitions. 

We now can proceed to derive an expression for the 
EAR, which equals to the rate in which energy goes from 
the system to the bath. Accordingly we have to calculate 
the rate of energy flow Q that is associated with the bath 
induced transitions: 



Q 



W^p n 



W^Pn 



A, 



(29) 



Note that in steady state this should equal the rate of 
work W, that is associated with the driving source induced 
transitions. We would like now to simplify the above ex- 
pression. The first step is to rewrite a\b\ — 0,2 &2 as the 
combination of (a\— a2)(&i+&2) and (ai+<X2)(&i— 62)- This 
gives 



Q 



E 



[(w<n-Wlt)Pn + «4(p«-Pn-l) 



(30) 



For the calculation of the differences in the above expres- 
sion we use Eq. (0) and Eq. (j2"5|) respectively, and further 
simplify by assuming Tb ^> A, leading to 



p n w 



PnW-, 







Tr 



-A 2 



T, 



(0) 



(31) 



In the above expression it was convenient to keep using 
the notation T^ ' even if the current is non-zero. 

The last term in the EAR expression Ea.([3~Tj) arises be- 
cause we have a non-trivial topology that can support non- 
zero current. We therefore name it 'topological term'. We 
see that this term depends on the ratio w^/(w^ + w^). 
Using the definition Eq. (fT^)) and the expression for the 
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Fig. 4: The energy absorption rate (EAR) versus the driving 
intensity for a Ring with N = 10 6 sites. The other parameters 
are the same as in the previous figure. The thin solid curve is 
the total EAR, while the dotted horizontal line is the expected 
saturation value. The decorated solid curve is the EAR differ- 
ence if the Ring is disconnected at one point, while the dashed 
curve is the topological term EqI. The current in the ring is 
also drawn (triangular markers). There is clearly a correlation 
between the excess EAR and the current. 



SMF Eq. fTo]) it takes the form Tb£qI- Consequently we 
can write the expression for the EAR schematically as in 
Ref. with an the additional topological term: 



Q 



Db 
T b 



Db 
T(°) 



T b £q I 



(32) 



The first term represents the heat flow due to a tem- 
perature gradient. The definitions of the diffusion coef- 
ficient Db and the induced temperature T^ are implied 
by a comparison with Eq. ()31l) . 

It is interesting to point out that Eq. (|3"2"|) connects be- 
tween the entropy production obtained under different lev- 
els of coarse graining [7]. The term SqI is a coarse-grained 
entropy production for a setup in which one can not distin- 
guish between transitions mediated by the thermal bath 
and noise. In contrast, Q/Tb is the entropy production 
rate when the two types of transitions are distinguishable. 

EAR vs Current.— We would like to inquire whether 
the EAR is correlated with the current. In a superficial 
glance one may be tempted to say that there is a linear 
correlation due to the topological term cx / in Eq. (|52"|) . 
However, one should realize that also the p n in Eq. ([3"TT) 
depend implicitly on the current. Hence one should not 
expect a strict linear relation, but, speculatively, a weaker 
correlation. To confirm this conjecture we calculate in 
FigQ]thc difference between the EAR of a connected Ring, 
and the EAR of the same ring after it had been discon- 
nected at one point. We indeed observe that there is an 
unambiguous correlation. 
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The driving induces non- vanishing current /, hence the 
EAR of a closed ring is larger than that of a linear chain. 
In the linear response regime one may expand the EAR 
in powers of e 2 . One should observe that the / related 



corrections are of order e 
pression: 



Q 



E 

n 

Db 



PnW^ 



Let us write the explicit ex- 



T B S I (33) 



[( 5n 6 2 )-(.g Iie 2 ) 2 +Var(.g)6 4 



(34) 



When comparing a closed ring to a disconnected ring, one 
should be aware that the last term, which we call 'topolog- 
ical term', should be either included or excluded respec- 
tively. Compared with a disconnected ring the topological 
term obviously increases the EAR, but from the expression 
above we see that the global dependence of the EAR on 
the driving intensity remains sub-linear. Namely, one real- 
izes that without the topological term the e 4 component in 
the EAR expression has the coefficient — g 2 , while with the 
topological term the net coefficient becomes —~g^ 2 , which 
implies larger EAR but still sub-linear. 
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Summary.— The study of transport in network sys- 
tems has numerous applications, notably in physical chem- 
istry, where the dynamics is commonly described by a rate 
equation. There is much interest in studying NESS cur- 
rents that are induced either by periodically varying sys- 
tem parameters |12j . or by stochastic driving. Assuming 
the latter, we have considered the NESS of a driven ring 
that is coupled to a bath, and found both the steady state 
current and the EAR. We have demonstrated how the 
ring-like topology and the sparsity lead to a glassy NESS 
with a non-trivial current dependence, exhibiting an in- 
teresting crossover from linear-response to saturation via 
an intermediate Sinai-type regime. 
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